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Abstract 



We develop and experimentally verify a theory of Type-II spontaneous 
parametric down-conversion (SPDC) in media with inhomogeneous distribu- 



tions of second-order nonlinearity. As a special case, we explore interference 

&\ 

effects from SPDC generated in a cascade of two bulk crystals separated by 
an air gap. The polarization quantum-interference pattern is found to vary 
strongly with the spacing between the two crystals. This is found to be a 
cooperative effect due to two mechanisms: the chromatic dispersion of the 
medium separating the crystals and spatiotemporal effects which arise from 
the inclusion of transverse wave vectors. These effects provide two concomi- 
tant avenues for controlling the quantum state generated in SPDC. We expect 
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these results to be of interest for the development of quantum technologies 

and the generation of SPDC in periodically varying nonlinear materials. 
PACS number(s): 42.50.Dv, 03.65.Bz, 42.65.Ky 
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I. INTRODUCTION 

Spontaneous parametric down-conversion (SPDC) [[[]] has now come into widespread use 
as a simple and robust source of entangled photon pairs. Uses for these pairs range from the 
examination of quantum mechanical foundations to applications in optical measure- 

ments [§], spectroscopy ||, imaging J7|, and quantum information fUJ. As such, there has 
been considerable interest in greater optimization and control of the exòtic two-photon states 
available from SPDC, particularly when pumped by ultrafast pulses |Ï0Hi2|. Additionally, 
much work has recently been focused on the use of cascaded nonlinear crystals [|Ï3H Ï8|1 to 
manipulate and improve the generation of the two-photon state. 

The photon pairs from Type-II SPDC are generated in a quantum state which can be 
entangled in frequency, wave vector, and polarization. In recent works |üfl , we have demon- 
strated the utility of a model which considers entanglement in these parameters concurrently. 
It was shown that quantum-interference patterns were altered predictably by controlling the 
range of transverse wave vectors selected by the optical system. In this paper, we extend this 
formalism and investigate interference from SPDC generated in media with inhomogeneous 
longitudinal distributions of nonlinearity. The state function of the photon pair generated 
in SPDC is completely characterized by three functions: the spectral profile of the pump, 
the longitudinal distribution of nonlinear susceptibility, and the dispersion in the generation 
medium. In principle, one could arbitrarily weight the spatiotemporal distribution of signal 
and idler modes by a judicious choice of these three functions. This, in return, introduces 
new avenues of control. To demonstrate this, we consider a rudimentary case of an inho- 
mogeneous medium, two bulk crystals separated by a linear medium such as an air gap. In 
this configuration, a host of interesting effects emerge, such as the modulation of interfer- 
ence visibility with crystal separation. This effect and others are theoretically predicted and 
experimentally verified in this paper. 
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The promise of a source whose degree of entanglement is controllable in freqüències 
and wave vectors by turning a single knob is clearly alluring for purposes of quantum- 
information processing. The results reported in this paper are also likely to be of use 
in guiding future developments in quantum-state synthesis involving multi-crystal configu- 
rations |T^-|TÜ[, in ultrafast-pumped parametric down-conversion [|I^,|n[], and periodically 
poled materials [pO,pT[ . 



II. THEORY 

Our theory considers a quantum state which can be concurrently entangled in polar- 
ization, frequency, and transverse wave vector, so as to be vàlid for an arbitrary optical 
system. As we shall see, the longitudinal distribution of nonlinearity provides a powerful 
means for controlling the structure of the two-photon quantum state generated in SPDC. As 
an important special case, we consider the simple example of SPDC generation in a cascade 
of two bulk crystals separated by a linear medium. We then describe the quantum inter- 
ference between the two photons of the SPDC pair as they propagate through an arbitrary 
linear optical system. This formalism allows the quantum interference to be analyzed in the 
absence of spectral filters and reduces to the conventionally established single-mode theory 
in the small-aperture limit, unless very thick crystals are used. 



A. State Generation in Inhomogeneous Media 

For the sake of simplicity, we consider media where effects from third- and higher-order 
susceptibilities are weak and can be neglected. By virtue of the relatively weak interaction 
in the nonlinear crystal, we consider the two-photon state generated within the confines of 
first-order time-dependent perturbation theory. The two-photon state at the output of the 
nonlinear medium is found in the interaction picture to be |T^] 

|^ {2) ) oc J dq dq e du du e $(q G , q e ; u , u e ) àl(q ,u ) àl(q e ,u e ) |0) , (1) 
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where the state function 

$(q ,q c ;^ ,^ e ) = È p (q + q c ; u Q + iü e ) J dz X (2) (z) e iA ^^> . (2) 

Here x^ 2 \ z ) is the distribution of second-order nonlinearity along the longitudinal axis, 
-Ép(q p ; u p ) is the complex-amplitude profile of the pump field, q., (j = p, o, e) is the transverse 
component of the wave vector kj in the medium, and A is the wave vector mismatch function 

A(q Q , q e ; u Q , u c ) = k p (uj + u c , q Q + q c ) - k d (u , q G ) - K e (u e , q e ) , (3) 

which depends on the dispersiveness of the medium. In this equation, the longitudinal 
projections Kj (j = p, o,e) are related to the indices (qj,Uj) via 

kj (ujj = yjtfiu^qj) - \cy\ 2 , (4) 

where uj v = u + oj c and q p = q Q + q c - Here the wavenumber kj = |k 3 -| = n[uj, 0(qj)] ujj/c, 
where c is the speed of light in vacuum, 9 is the angle between k p and the optical axis of the 
nonlinear crystal, and n(cüj, 9) is the index of refraction in the nonlinear medium. Note that 
the symbol n(o>, 9) represents the extraordinary refractive index n e (u>, 9) when calculating k 
for extraordinary waves, and the ordinary refractive index n (uj) for ordinary waves. 

Note from Eq. @ that the state function is completely characterized by the spectral pro- 
file of the pump, the longitudinal distribution of nonlinear susceptibility, and the dispersion 
in the generation medium. All three of these parameters may be controlled experimentally, 
and all three present avenues for controlling the structure of the two-photon quantum state. 

For a medium with an inhomogeneous distribution of nonlinear susceptibility along the 
longitudinal axis, it is convenient to define 

X {2 \z)= Jd(x (2 \Oe-* z , (5) 

where x (C) is the inverse Fourier transform of X ( z )- Substitution into Eq. (^|) then gives 

<£>(q ,q e ;^ ,^ e ) = É p (q + q e ; u Q + u c ) x (2) [A(q Q , q e ; u Q , u e )] (6) 
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for a uniformly dispersive medium. For example, a single bulk crystal of thickness L and 
constant nonlinearity Xo has a nonlinear susceptibility profile x^ 2 \ z ) — Xo Tec ^[-L,o]( z ) where 
rectr_i i Qi(^) = 1 if —L < z < and zero otherwise. In this case, the inverse Fourier transform 
of the nonlinearity profile becomes 

X (2) (A) = Xo^sinc(^)e-^. (7) 

For a monochromatic plane wave pump with a central frequency u p , E p (q p ;uj p ) in Eq. @ 
is proportional to <5(q Q + q e )S(uj + u e — u p ) and the state function $ for SPDC reduces to 



X (2) 



A(q, — q; lo, — uj) . Figure l(a) shows the absolute square of the state function in 



p 

■2 



Eq. (0), the familiar sine (LA/2) distribution of SPDC from a single bulk crystal. 

1. Periòdic Nonlinearity 

We now consider a medium of thickness L with a periòdic distribution of nonlinear 
susceptibility x^ 2 \ z ) — X ( z + A) within the medium. Such materials are widely used in 
classical nonlinear òptics |2(J and have recently been employed for generation of SPDC [21 



We may write 

X (2) (z) = Xo9( z ) rect[_L j0 ] ( z ) , (8) 
where g(z) can be expressed in the Fourier series 

oo 

g(z)= Yl Gme iK - z (9) 

m=— oo 

with K m = 27rm/A. The Fourier Transform of Eq. (§) is then given by 



X^(A) = Xo L £ Gm 



oo 



sine 



^(A + K m ) 



For example, let us consider the case of a sinusoidal distribution of nonlinear susceptibility 
with period A, for which 

X {2 \ z ) = Xocos (j^A nxb[-L,o](z) , (11) 



which yields 



X (2) (A) = XoL {sine [f (a + f )] e^ï) + sin c [f (a - |)] e^^)} . (12) 

In this case, we obtain phase-matching conditions similar to the first-order quasi-phase 
matching (QPM) observed in periodically poled nonlinear crystals. The extra component 
±27r/A above is analogous to the grating vector in first-order QPM. Figure l(b) shows the 
absolute square of the state function of the down-converted light obtained from a single 
crystal with the nonlinearity profile given in Eq. (ITT]). 



2. Cascaded Bulk Crystals Separated by Linear Media 

A simple example of a medium with an inhomogeneous distribution of nonlinearity is a 
cascade of múltiple bulk crystals separated by linear dielectrics. Consider, for example, a 
cascade of A" bulk crystals separated by A^ — 1 linear media. Let each nonlinear crystal j 
have thickness Lj, constant nonlinearity \0ji an d separation distance dj from the previous 



crystal. The overall nonlinear susceptibility of this system is then given by 

N 

,0] 



X {2 \z) = e jX0j rect [ _ L . jl 



N 

k=j+i 



(13) 



3=1 

where the terms of the summation inside the rect function are taken to be zero if k > N. Here 
e = ±1 represents the sign of the quadratic susceptibility, which depends on the orientation 
of the optical axis of the j-th crystal. Note that in this equation, the z = point is placed at 
the output plane of the last crystal. In such a configuration, the function x i n Eq. (^) 
becomes 

X (2) (A) = f; e jXOj L sinc (L·l^£\ ^££1^ (14) 

j= i 3 V 2 / 

where the wave vector mismatch function Aj is independent of z and A'. As seen from 



Eq. (^|), A' depends on the dispersiveness of the linear medium [Ï4| . 

We now consider the particular case of a cascade of two bulk crystals of the same ma- 
terial [See Fig. 2(a)] separated by a linear but dispersive medium. The explicit form of the 
nonlinearity is given by 
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X {2 \z) = Xorect[_ Llj0 ](2: + d + L 2 ) + e Xo rect[_ L2i0 ](2;) , (15) 

where e = +1 if the optical axes of the two crystals are parallel and e = — 1 if the optical 
axes are antiparallel. For such a configurat ion 

X (2) (A) = Xo {^isinc (í^) e-^e^AW) + ^ e -i^j (16) 

The absolute square of the state function of SPDC in this configuration is given in Fig. 2(a), 
where the envelope is governed solely by the dispersion in the nonlinear crystals. The period 
of the modulation inside this envelope is determined primarily by the dispersion in the linear 
medium between the crystals, while the amplitude of this modulation is determined by the 
ratio of the crystal thicknesses. Figure 2(b) shows the absolute square of the state function 
in the special case of two bulk crystals of the same material with the same thickness. In this 
condition the amplitude of the modulation inside the envelope is maximized. 

B. Two-Photon Amplitude and Fourth-Order Correlation 

We now consider the propagation of the down-converted light through an arbitrary linear 
optical system to a pair of detectors, as illustrated in Fig. 3. The joint probability amplitude 
of detecting the photon pair at the space-time coordinates (x A ,£ A ) and (x b ,íb) is given by 

A(x a ,í a ;x b ,íb) = (0|4 +) (x a ,ía)4 +) (xb,íb)|^ (2) ) (17) 

where and are the positive-frequency components of the elèctric fields at points A 
and B. The explicit forms of the quantum fields present at the detection locations are given 
by 

£ A +) (x A ,í A ) = E fdqdw e^ ÍA W A j(x A ,q;w) %(q,^) , 

j=e,o 

4 +) (x B , Í B ) = E / d ^ duü e "" <B W Bj -(x B , q; u) àj(q, oo) , (18) 

j=e,o 

where the transfer function Hij (i = A, B and j = e, o) describes the propagation of a mode 
(q, cu) through the optical system from the output plane of the nonlinear medium to the 
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detection plane. Substitution of Eqs. (HD and (|Ï8|) into Eq. (|Ï7|) yields a general form for the 
two-photon detection probability amplitude, 

^4(XA, Ía; Xb, Íb) = ^4ao,Bc(xa, Ía! Xb, Íb) + ^Bo,Ae(xA, Ía! Xb, íb) (19) 

where the probability amplitude A AoBe for finding the signal photon in arm A and the idler 
photon in arm B is defined as 

^Ao,bc(xa,ía;x b ,íb) = / dq dq e düj du e $(q Q , q e ; oj 0} u e ) e -^otA+^t B ) 

x7íao(x a , q ; w ) ^bc(x a , q c ; w e ) (20) 

and y4Bo,Ae(xA, ía; x b, íb) is obtained by exchanging the indices A <-» B. 

The joint probability density for detection of the signal and idler photons at space-time 
points (xa, ía) an d ( x b, íb) is given by the fourth-order correlation function, is given by the 
absolute square of Eq. ([Ï9| ) : 

, Ía;X B ,Íb) = |^Ao,Be(xA, ÍA! X B , Íb)| + |^Bo,Ac(xa, Ía! Xb, Íb) | 

+2 3fíe[v4X 0iB e( x A, ía; x b , í b ) v4 Bo ,a c (xa, *a; x b , íb)]- (21) 

With current technology, quantum interferometry is performed using slow detectors that 
cannot resolve signals on the characteristic time scale of down-conversion (the inverse of 
down-conversion bandwidth), which is typically less than 1 ps. In addition, the detectors 
used in our experiments have a large active area compared to the width of the SPDC beams 
at the detection planes. Under these conditions, the coincidence count rate R is readily 
expressed in terms of the two-photon detection probability amplitude A by integrating the 
fourth-order correlation function G^(xa, ía; xb, íb) over all space and time, 

R = J dt A dt B dx A dx B ^(xa^aïxb^b)! 2 - (22) 

This expression for the count rate can be separated into two terms as 

R = Ro + -Rint , (23) 
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where the baseline term is 



Rq = J du' du dq Q dq_ e dq^dq^ $(q D , q e ; U, u') $*(qó, qó; W, tu') 

x [5 A B(qo,qc,qó,qc;^,^') + 5 BA(qo,qc,qó,qé;^^')] ( 24 ) 

and the interference term is 

fíiüt = 2Ke J du' du dq Q dq e dqó<iqó $(q Q , q e ; u/, tu) $*(qó, qó; u/) SAB(q c , q , qó, qó; d/) . 



In Eqs. (0) and (25) the state function $ weights the signal and idler modes in the process 
of generation, while the function Sab weights these modes in the process of propagation 
through the optical system. Explicitly, 

5 A B(q c , q , qó, qó; <*>') = (^a ( x a, qó; w) 7í A o(x A , q G ; w)) XA (ft Be ( x B, qó; w') ^b c (x b , q e ; w')> 

(26) 

where (-) Xj indicates integration over the total detector area. 

Note from Eqs. (|T9|) and fl20|) that the two-photon detection probability amplitude is 
completely specified by Tïij (i = A, B and j = o,e), <£>(q , q e ; u; , u; e ), and the physical 
location of detectors A and B. As we have seen in Eq. (^|), we may control the structure 
of the state function $ by a judicious choice of the pump spectral profile, the longitudinal 
distribution of nonlinearity, and the dispersion in the crystal. We may further control the 
two-photon detection amplitude, and hence the quantum-interference pattern, through the 
choice of the optical system. Note that states with different state functions can lead to the 
same quantum-interference pattern through an appropriate design of the optical system. 

In the experiment ally relevant case of a monochromatic plane wave pump field, Eqs. (p4|) 
and (125D become 

Ro = J du J dq dq' x (2) * [A(q', -q'; u, u° p - u)} x (2) [A(q, -q; u, u° y - u)} 

x [Sab (q, q'; w) + Sba (q, q'; w) (27) 

and 
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Bint = me I dqdq' x (2) * [A(q', -q'; u° p - cu, cu)] x (2) [A(q, -q; u, Lü° p - cu)] S AB (q,-q» (28) 

where we use the shorthand 

^AB(q, q'; üj) = <Sab (q? -q, q\ -q 7 ; w, wj - w). (29) 

Thus we see that for a monochromatic plane wave pump, the quant um-interference pattern 
is critically dependent on the form of x^ 2 \A), which we are free to choose as a design 
parameter ||20|| . In principle, the only limitation on the class of amplitudes A which we are 
able to prepare with this method is the restriction that the optical system is linear. 

C. Quantum Interference with a Cascaded Pair of Bulk Crystals 

We now apply the above formalism to the case of two cascaded bulk crystals separated by 
a dispersive but linear dielectric medium such as an air gap. For simplicity, we again consider 
the medium to be pumped by a monochromatic plane wave. Owing to the structure of the 
nonlinearity for this particular case [Eq. fllg)], the overall two-photon detection probability 
amplitude is the sum of the two amplitudes associated with each single crystal |fL4|,|Ï8l| . Each 
of the amplitudes in Eq. (|Ï9"D can then be written as 

^Ao,Bc(xa, X B , Íb) = ^Ao,Bc( x A ; t A ] X B , £ B ) + 6 ^. A o,Be( X A> ^A! X B , Í B ) (30) 

where e = ±1 as in Eq. (|Ï5|) , and a similar expression for A-Q ,Ae is obtained by exchanging 
the indices A <-» B. ^AoBe f° r r = (^2) is the probability amplitude of finding the o- 
polarized photon generated in the r-th crystal in arm A and the e-polarized photon from 
the r-th crystal in arm B. From Eq. (|20|), 

^Ao,Be( x A,x B ; í) = J dudqe~ lut x^\cí, v) 7í A o(x A ; q, v) H Be (x B ; -q, -v) (31) 

where the angular frequency v = uj — is the deviation from the central frequency 

t = t a — t B is the time difference between detection events, and (q, v) is the inverse 

(r) 

Fourier transform of the nonlinearity profile of the r-th crystal. A Bc , Ae is likewise obtained 
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by a suitable exchange of the indices. Note that we have omitted an overall phase factor 
exp[— zcijp(ÍA + Íb)/2] which appears outside the integral in Eq. (|3Ï|), since in experimental 
practice we are interested only in the absolute square |A(xa, ía; x b, íb) | 2 , so this factor does 
not introduce any relative phase between the terms of Eq. (|T9|). 

As we are considering bulk crystals, the nonlinearity profile of each crystal r is uniform, 
and thus 

Xl 2) (q, v) = XoQr (q, u)-Jdz rect hLri0] (z) e iA ^> , (32) 

where 



Q r ( q , v ) = e -i[dA'(q^)+i 2 A 2 ( qiI ,)] 5 r ,i (33) 

is the transfer function for propagation of the signal (o-polarized) and idler (e-polarized) 
fields generated in the first crystal through the linear dispersive medium of thickness d and 
thence through the second crystal of thickness L 2 . Alternatively, Q r (q, v) may be thought 
of as the phase accumulated in the shift of the rect function for the first crystal by a distance 
— (L 2 + d). As given in Eq. (|3D, A r is the wave vector mismatch function due to dispersion 
in the r-th crystal and A' is the wave vector mismatch function due to dispersion in the 
linear medium. The symbol 5 r i represents the Kronecker delta where 5n = 1 and 82 1 = 0. 



1. Coíncidence Detection 

Taking the absolute square of Eq. (|30| ) gives interference between the probability am- 
plitudes of finding a pair generated in the first crystal and finding a pair generated in the 
second crystal. Indeed, substitution of Eqs. (|Ï9| ) and ( |30| ) into Eq. (|22| ) gives the coincidence 
count rate as a sum of three contributions 

R = + R^ + R^ (34) 

where the first two terms are the coincidence-count rates for single-crystal SPDC, and the 
last term arises as interference between the two single-crystal amplitudes. Recalling Eq. (PB|). 
each term in Eq. ( ^4|) can in turn be broken down into baseline and interference terms 
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R (h) = R (h) + r: 



(35) 



where h = 1,2, 12. 

We now use this theory to predict quantum-interference patterns in the interferometer 
shown in Fig. 4. In this system, the transfer function 7i%j is separable into diffraction- 
dependent and -independent factors as 



where the polarization-independent components are grouped into Hi and the remainder are 
grouped into %j é Ki ^ v > lr . In this case, %j = (e^ ■ e^) is the projection of the unit photon 
polarization vector e,- (j = o, e) onto the axis of the polarization analyzer in front of detector 
i = (A, B), and the exponential factor is the transfer function of the delay line. 

The delay line, which is often treated as a simple phase shift, is a dispersive optical 
element which may alter the spatial and/or the spectral profile of the two-photon probability 
amplitude. Experimentally, the delay line consists of a birefringent quartz plate of variable 
thickness, modelled by the propagation function e tK3 ^^ lT , where the longitudinal projections 
Kj of the signal and idler wave vectors are defined in Eq. and l T is the thickness of the 
birefringent plate, which induces a relative optical-path delay r. This propagation function, 
while here used to describe the delay line, is technically a vàlid transfer function for any 
nonabsorbing dispersive optical element. 

Given this particular optical system, then, the single-crystal coincidence-count rates are 
given by 



TCij^Cj-, q, v) = Hi( 



(36) 




X ^Ao.Bc X 



f*(q', v) xl 2) (q, v) + /4o,Ae } *(-q\ -v) X?\~ci, -v) 



(37) 



and 




X /iAo.Be /^Bo,Ae Xrtà, ~ V ) Xr fa» V ) + Xr , v) Xr ("Q, ~ v ) 



(38) 
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where r = 1,2 is the crystal index, \ii^i m = %j ■ % m = A, B and j, m = e,o) is the 
projection of the polarization of the i-th photon onto the j-th basis polarization and the 
polarization of the Z-th photon onto the m-th basis polarization, and the phase function 



T] T (q, v) = - [k (u, q) + K e (-u, -q)] l T 



(39) 



depends on the dispersion introduced by the delay line. The integral over the detection 
planes 

JF AB (q, ±q' , v) = (H A (x A , q; v) H* A (x A , ±q'; v) ) (H B (x B , -q; -v) (x B , =Fq'; -v) ) (40) 



is an analog of the function iSab for the polarization-independent elements of the system 
only 

Meanwhile, the coincidence-count rates which arise collectively between the contributions 
from the two crystals are given by 

= e f dudqdq' F AB {q,q',v) e-^^-nrW *)] 

■oPW»' r A ,~( 2 ) ( n ,A J_ „2 v( 2 W_«' _,A ,~(2) 



r=l,2 



(41) 



and 



^int = e / du dqdq' Fab(ci, -q', ^) e" 



-i[r? r (q,^)-í? r (q',-i/)] 



/^Ao,Be A'Bo.Ae 



r,(2) 



;q^) + xl 2) *(-q'^)# r (-q,-^)l ■ (42) 



r=l,2 



By combining Eqs. (0), fl38|) , (|4l"|), and (42), the overall coincidence-count rate -R(t) can 
be organized into a general form 



R(t)=R [1+v po1 V{t)\ 



(43) 



where Rq is the baseline coincidence-count rate and the overall projection of both photon 
polarizations onto the basis of the polarization analyzers is given by the factor 

/^Ao,Be /^Bo,Ae 



^pol 2" 2 2 

^Ao.Be + /^BojAe 



(44) 



Observe that the r-dependence of the quantum interference pattern is then contained solely 
in the visibility function V(r). 
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2. Cascade of Two Identical Crystals 



In the experiments presented in this paper, the two crystals are of the same material and 
have equal thicknesses L. The apertures are symmetric for both transverse directions, and 
the analyzers are set 45° from the optical axis, so v po \ = — 1. No spectral filters are used. 
For these conditions, the explicit form of the visibility function in Eq. fl4*3|) becomes 

V(r) = — ^ / dz U L (z) Yl L ( 2 4--2L-z) G {1) 1 ' 7 



l + p 2 J w \D ) \V LD 

+-r^ í dz U L (z) n L (% - z) ' : T 



l + p 2 J \D J W LD 



[ L 2'LD' 



(45) 



where is defined in Appendix ^|, p = (di + d)/di, ILl(z) is the unit rect function from 
[0, L], and D = u^ 1 —u^ 1 is the dispersion coefficient of the nonlinear medium. It is through 
the £/-functions that spatiotemporal effects enter the quantum interference pattern. Details 
on the derivation of Eq. ([Ï5|) can be found in Appendix [Al. 

The collective interference term in Eq. (|45|) shows interesting behavior in certain limits 
of crystal separation. If the optical axes of the two crystals are parallel, the coincidence- 
count rate reduces to that from a single crystal of thickness 2L as d — > 0. Further, it 
reduces to that from a single crystal of thickness L as d — >• oo. We also note the absence of 
any shoulder modulation with r, an important indication of the purity of the polarization 
Bell-state which is formed in postselection. In the case where the two crystals have equal 
thickness, the strongest interference occurs at delay r = LD. The visibility at this point is 
given by 

V(LD) = 2e-^— 2 i' d( Xe [G {12) (C; 1)] (46) 
1 + p À Jo 1 J 

where ( = z/L is a convenient dimensionless variable. 
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3. Small-Aperture Approximation 



In the limit of very small apertures, no transverse wave vectors are allowed to propa- 
gate through the interferometer, and in the case of sufficiently thin bulk crystals and small 
separation distances, the quantum interference is effectively described by the convent ionally 
used single-mode theory [|Ï%P!2|| . 

Figure 5 is a sketch illustrating how quantum interference arises in the interferometer of 
Fig. 4 assuming two identical crystals of thickness L\ = Li = L, dispersion coeficient D, 
paral·lel optical axes, and both polarization analyzers set to 45°. In the limit of sufficiently 
small apertures, we may apply the conventional single-mode theory |22| and write the third 
term in Eq. (|2Ï| ) as the product of two probability amplitudes A(t^ — íb) and A*(— ía + íb), 
which slide back and forth across the ía — axis as the relative optical-path delay r is varied. 
In the diagonal portion of the illustration, these amplitudes are depicted by two grey-and- 
white rectangles. When the delay is set such that the two rectangles overlap, interference 
can be seen. 

Within each rectangle, the white box represents the probability amplitude for detecting 
a photon pair produced in the first crystal, while the grey box represents the detection 
amplitude for a pair produced in the second. As the delay is set in the region < r < 
LD/2, the grey boxes overlap (shown as black) but the white boxes do not. In this regime, 
interference typical of a single crystal of thickness 2L is observed. As the delay is increased 
into the region LD/2 < r < 3LD/2, the photon pairs produced in the first crystal become 
indistinguishable with the photon pairs produced in the second crystal at the detection 
planes. As such, the probability amplitudes of detecting photon pairs produced in each each 
crystal exhibit collective interference. This is seen pictorially by the overlap of the grey 
boxes with the white boxes. Note that the two rectangles overlap completely at r = LD, 
the center of the region of interference. The interference visibility in this region depends 
on the phase shift between these two amplitudes, which is in turn dependent on spatial 
effects and the dispersion in the linear medium separating the crystals. As the distance d 
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between crystals is changed, the visibility in this region modulates sinusoidally between ±1. 
As the delay is increased still further into the region 3LD/2 < t < 2LD, the white boxes 
overlap (shown as light grey) but the grey ones do not, and we again return to the regime 
of single-crystal interference. In this way we can trace out the interference dip predicted by 
Eq. P5|). This is illustrated in the inset at the lower right while representative experimental 
data is shown in the upper left inset. Figure 6 is similar to Fig. 5, but for the case of unequal 
crystal thicknesses L\ ^ L 2 . 

In this small-aperture approximation, the dispersion-related ^-functions in Eq. ( f45[) for 
the cases of parallel (p) and antiparallel (a) optical axes at t/LD = 1 are found in Appendix 
0to be 

"(C-i) 2 (C + i) 2 ' 



^ 12) (C; 1) = expHA'd) exp Í^J^l 



and 



ÇW(Ç- 1) = exp{-tA'd) exp <| -ï^ 



L·lMLl 2 



d\ di + d 



(C-i) 2 (C-i) 2 



(47) 



(48) 



di di + d 

where k p = 2ir/\ p and |M| = ln[n e (c<jp/2, 9 P )] is the spatial walk-off vector [23 . 

When a plane wave pump is normally incident on thin crystals and the apertures are 
sufficiently small, the strictly collinear signal and idler beams are selected by the optical 
system, and no spatial effects due to transverse wave vectors can be observed. In this 
limit the dominant contribution to the phase between probability amplitudes arises from 
dispersion in the linear medium between the crystals [see Eqs. (f47|) and (f48|)1. For degenerate 



SPDC this phase term in air is [24 



<p disp (d) = A'd = k p [n(X p ) - n(2A p )] d ~ tt 0.059 d[mm] (49) 

where we have taken À p = 351.1 nm in keeping with our experiments. 

However, when long crystals are used, or the separation distance between the two crystals 
becomes large, some spatial effects due to transverse wave vectors become observable. The 
visibility at the center of the region of interference tau = LD [Eq. fl4"E|)] for the case of 
parallel optical axes is given by 



V P (LD) = 2 



dl + (di + d) 2 Jo 



j dÇ cos < 


í K 


ML 


2 " 


2(di + d) 



+ (pdi sp {d) 



(50) 

while for the case of antiparallel optical axes, it is given by 

' (51) 



k p \ML\ 2 d 2 

(C - 1) + 0disp(^) 



2(d 1 + d) d x 

Furthermore, for small separation (d — > 0) between the two crystals with parallel optical 
axes, Eq. ( f46[) becomes 

V P (LD) = sine (J^-\ML\ 2 ^j (52) 

as opposed to the unity value predicted by the single-mode theory. Note that this visibility 
is identical to that of a single crystal of thickness 2L as was shown in previous works 



Ï2| . Applying parameter vàlues typical of our experiments, we find that if d± = 1 m, 
À p = 351.1 nm and |M| = 0.07, then V p ~ sinc(0.044 í[mmf). For these conditions, 
therefore, an observable deviation of V p from unity can be realized only for sufficiently large 
crystal thicknesses. 

In the case of antiparallel optical axes, the visibility at the center of the region of inter- 
ference is 

V a (LD) = -l, (53) 

which means that the coincidence-count rate is twice as high at the center as it is on the 
shoulders, assuming the analyzers are set so that v po \ = — 1 [see Eq. ([43])]. In this case, the 
minus sign in Eq. (|53D arises from a sign difference between the quadratic susceptibilities 
of the two crystals, and hence a sign difference between the spatial walkoff vectors Mi and 
M 2 . In effect, the spatial walkoff in one crystal compensates for the spatial walkoff in the 
other, and the spatial effects due to transverse wave vectors are cancelled out. Exploitation 
of this effect is common in the design of optical parametric oscillators. 
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4- Multi- Parameter Formalism: Spatial Effects 



The quantum state generated from SPDC, which is concurrently entangled in frequency 
and transverse wave vector, leads to transverse spatial effects that can be observed in quan- 
tum interference. As the aperture diameters are increased or the apertures are brought closer 
to the output plane of the nonlinear medium, a greater range of wave vectors is allowed to 
propagate through to the detectors. These transverse wave vectors introduce distinguisha- 
bility between the signal and idler photons, thus reducing the visibility of the observed 
quantum interference. This is analogous to the temporal distinguishability introduced by 
the use of a femtosecond-pulsed pump. 

Equation fl45|) is vàlid for any linear optical system. However, to enable swift evaluation of 
the integrals in this equation we approximate the circular apertures used in the experiments 
by "soft" Gaussian apertures of (1/e) widths and r B . A sharp circular aperture, of the 
type used in experiments, has a diffraction pattern described by a first-order Bessel function, 
whereas a Gaussian aperture, of the type used in the numerical simulations, has a Gaussian 
diffraction pattern. Despite this fundamental difference, it is a fair approximation if the 
width r of the Gaussian is selected to roughly fit the width b of the Bessel function. In our 
calculations, this is done by choosing r = b/2 \/2. This approximation offers an indispensable 
advantage, as it allows us to evaluate the ^-functions analytically (see Appendix [B]) and 
thereby reduce the demand for numerical integration in making theoretical predictions. 
Under this approximation, an expression for the visibility function at (r = LD) for parallel 
orientations of the optical axes is given by [see Appendix |TJ and Eq. (p6j)] 

di(di + d) 1 



VJLD) = 2 



where 



d\ + (d 1 + df v^+^F 

x f 1 dC e'^^^k) cos [cC 2 + X>C + £-0 7 (d) + 0disp(d)l (54) 
Jo 1 J 



k v d 

7 = h j I a ' 07 ( d ) = - arctan( 7 ) , (55) 



4c?! di + d 
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and r 2 = [r\ + r|)/2. Here 

B = 2 
C 
V 
£ 



k p \M\Lr 
4(di + (i) 



fc p |ML| 2 d 



1 / _dj 
T+^[ + 2d lj 
1 



8(0?! + d) dil + 7 2 



A; P |ML| 2 
~2(di + d) 
A; P |ML| 2 



d 



1 



d 



2di/ 1 +7 2 

20?! + cA 1 
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(56) 
(57) 
(58) 
(59) 



^(di + d) \dt * ' d / 1 + 7 2 ' 
Figure 7 presents a plot of the modulus of V P (LD) given in Eq. ([54]) as a function of the 
crystal separation d for aperture diameters b = 2.5 mm, 4.0 mm, and 5.0 mm. Sòlid, shaded, 
and open circles denote the màxima of visibility and correspond to 2.5-mm, 4.0-mm, and 5.0- 
mm aperture diameters, respectively. For this plot, L = 0.5 mm, A p = 2ir/k p = 351.1 nm, 
and \M\ = 0.07 in keeping with our experiments. 

For antiparallel axes, the visibility function at r = LD is 

di(di + d) 1 



V a (LD) 



J d 2 + (d! + d) 2 vTT^ 



X 



dC e" 



-S(l-Cf 



cos 



^(l-C) 2 -0 7 (d) + disp (d) 



(60) 



For the case of two crystals in contact (d = 0) we have B = 2(A; p |M|Lr/4di) 2 , C = 0, 



V = -A; p |ML| 2 /2d 1 , and £ = 0. From Eq. 
t = LD for parallel optical axes is 



we now observe that the visibility at 



V P (LD) = exp 
while for antiparallel axes 



(k p \M\Lr 
~ 2 V 4dx . 



sine í -?-|MIr| 



(61) 



V a (LD) = -V2^ 



di 



k p \M\Lr 



erf 



A; p |M|Lr 
, 2^f2d x , 



(62) 



Note that these visibilities depend on the aperture diameter, as well as the crystal thickness 
(as was seen in the previous section), and are once more markedly different from V a = V p = 1 
as predicted by the single-mode theory. Equations ( |6"TD and ( |6"2"D are plotted in Fig. 8 as a 
function of Gaussian aperture width r for three different crystal thicknesses L = 0.5 mm, 
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1.5 mm, and 5.0 mm. As predicted, the plot of the visibility obtained with a parallel 
orientation of the crystal axes (sòlid) reduces faster than the visibility for an antiparallel 
orientation (dashed). 

III. EXPERIMENT 

A. Experimental Arrangement 

The experimental arrangement is illustrated in Fig. 9. A 200-mW cw Ar + -ion làser 
operated at 351.1 nm served as the pump. This highly monochromatic làser beam was 
passed sequentially through a cascaded pair of /3-barium-borate (BBO) crystals each with 
thickness 0.5 mm. The thickness of the air gap between the crystals was varied from d = 
2 mm to 100 mm. The crystals were aligned to produce pairs of orthogonally polarized 
photons in degenerate collinear type-II spontaneous parametric down-conversion (uf = uf = 
Up/2, where uf, uf, and u® represent the central freqüències of the signal, idler, and pump 
fields, respectively). The làser power was sufficiently low to ensure, with high probability, 
that at most one photon pair was generated at a given time. The high visibility obtained 
from separate single- crystal quantum-interference experiments confirmed the validity of this 
assumption. A dichroic mirror, which transmits the 702-nm signal and idler beams while 
reflecting the 351-nm pump, was placed after the two crystals to remove residual pump làser 
beam. 

The relative optical-path delay r was introduced using a z-cut birefringent crystalline- 
quartz plate of variable thickness. The range of transverse wave vectors for the down- 
converted light was selected by circular apertures of diameters 2.5 mm to 5.0 mm. These 
apertures were positioned at 750 mm from the output plane of the second crystal. The single 
aperture used in the experimental setup is formally equivalent to the use of two apertures of 
identical diameters in the interferometer of Fig. 4, as used for the theoretical discussions of 
previous sections. The beams of down-converted light were then directed to a nonpolarizing 
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beam splitter, and thence to the two arms of a polarization intensity interferometer. Each 
arm of the interferometer comprised a Glan-Thompson polarization analyzer set at 45° 
with respect to the horizontal axis in the laboratory frame, establishing the basis for the 
polarization measurements. This basis was selected so as to permit observation of the 
quantum-interference pattern as a function of the optical-path delay r. Finally, a convex Iens 
(not shown in Fig. 9) was used to reduce the beam size to be less than the area of the detector, 
an actively quenched Peltier-cooled photon-counting avalanche photodiode. No spectral 
filters were used in any of the experiments. Coincidence detection was performed using a 
3-nsec integration window and corrections for accidental coincidences were not necessary. 

B. Experimental Results 

First, we report quantum interference from crystals oriented with paral·lel and antiparallel 
optical axes in the small-aperture approximation. We demonstrate that the visibility at 
t = LD varies sinusoidally with crystal separation d. Second, we investigate spatial effects 
on the visibility at r = LD arising from the acceptance of a broader range of transverse 
wave vectors as the aperture is opened. 

1. Parallel and Antiparallel Optical Axes 

In this section we investigate the quantum interference patterns from SPDC in two 
identical crystals with parallel and antiparallel optical axes. The details of the experimental 
setup can be found in the previous section. The signal and idler fields were selected by a 
2.5-mm circular aperture positioned 750 mm after the second crystal. 

In the lower portion of Fig. 10, we plot the visibility function at the center of the inter- 
ference region, V(LD) as defined in Eq. (f4~6|), as a function of the separation d between the 
two 0.5-mm-thick BBO crystals. Each data point on the graph was obtained by calculating 
the ratio (Rq — Rmt)/Ro where Rç, and i?i nt are the coincidence-count rates on the shoulders 
(r < and r > 2LD) and at the center of the region of interference (r = LD), respectively. 
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The shaded circles plot the visibility at r = LD as a function of crystal separation d when 
the optical axes of the two crystals are parallel. The open circles correspond to the case of 
antiparallel optical axes. As anticipated by Eqs. (|54]) and (|60"D, the modulation in visibility 
with crystal separation is sinusoidal. The 7r-phase shift between these two cases arises from 
the change in sign of x as the relative orientation of the optical axes is inverted. The 
theoretical curves (sòlid) are plots of Eqs. (|4]) and (|60|) assuming a Gaussian aperture of 
width r (defined as the diameter b of the circular apertures divided by 2 a/2, the width at 
1/e), set to satisfy 2\^2r = 2.5 mm. The agreement of the theoretical curves with the 
experimental data demonstrates the validity of the Gaussian-aperture approximation. The 
remaining experimental parameters for this plot are d\ = 750 mm, À p = 2n/k v = 351.1 nm, 
and \M\ = 0.07. 

In the top left inset of Fig. 10, we show a representative example of quantum-interference 
patterns in which the visibility function V P (LD) for crystals with parallel optical axes is 
either minimum or maximum. In the pattern exhibiting a peak (minimum), the crystals 
were 17.5 mm apart, while in the pattern exhibiting a dip (maximum) the crystals were 
37.5 mm apart. Each data point presented on the visibility graph (lower inset) is extracted 
from such interference patterns. In the figure at the top right, where the two crystals are 
oriented with antiparallel optical axes, the opposite is true: the pattern exhibiting a dip 
(maximum) corresponds to a 17.5-mm crystal separation, while the pattern exhibiting a 



peak (minimum) corresponds to a 37.5-mm separation. Note from Eq. (f44|) that when both 
polarization analyzers are set to 45°, the polarization projection factor t> po i = —1. Thus, 
from Eq. (^), positive vàlues of V(LD) correspond to a quantum interference dip while 
negative vàlues correspond to a quantum interference peak. Orienting the crystals with 
antiparallel optical axes produces a sign difference between the second-order nonlinearities 
of the two crystals, so this condition is again reversed. 
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2. Spatíotemporal Effects 



In previous work with single-crystal SPDC [12|, we found that the acceptance of trans- 
verse wave vectors had substantial effects on quantum-interference patterns. To investigate 
spatíotemporal effects in dual-crystal SPDC, we carried out identical experiments to those 
detailed in the previous section, except that the 2.5-mm circular aperture was replaced with 
4.0-mm and 5.0-mm circular apertures. The apertures remained positioned 750 mm from 
the output plane of the second crystal. 

As with the single-crystal experiments, we found that increasing the diameter of the 
aperture lowers the visibility of the quantum interference patterns. This effect can be ex- 
plained by considering the two independent mechanisms that introduce phases between the 
two-photon probability amplitudes for each crystal: 1) dispersion in the air separating the 
two crystals, which results in the phase 0dis P [Eq. (|49|) 1, and 2) the angular spread of down- 
converted light which results in the phase 7 [Eq. (|55|)1. Whereas 0dis P is dependent on 
the crystal separation d, but independent of r and di, 7 depends on all three parameters. 
However, when d/di 3> 1, 7 ~ — arctan(/c p r 2 /4<ii). In this limit, 7 is dependent on r and 
di, but independent of the crystal separation d. As such, 7 dominates when the aperture 
is sufficiently large, while 0di sp dominates when the aperture is sufficiently small. 

Fig. 11 displays plots of the visibility function at r = LD as a function of the crystal 
separation d for 2.5-mm (top), 4.0-mm (middle), and 5.0-mm (bottom) aperture diameters. 
The two crystals are oriented with parallel axes. Note that for a fixed value of d, there 
is a reduction of visibility for increased aperture diameter. A reduction of visibility also 
occurs as the distance between the crystals is increased. Moreover, the period of oscillation 
contracts slightly as the aperture diameter is increased. 
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IV. CONCLUSION 



We have developed a theory of Type-II spontaneous parametric down-conversion (SPDC) 
in media with inhomogeneous distributions of nonlinearity. The down-converted light can 
be concurrently entangled in frequency, wave vector, and polarization. We have shown that 
the state function of the down-converted light can be controlled by design of the nonlinearity 
profile in the crystal, as well as the spatial and spectral profile of the pump field. As one 
rudimentary design, we have considered the case of two nonlinear crystals separated by a 
linear dispersive medium. We anticipate that even greater control can be identified and 
exploited when other distributions of nonlinearity are employed. 

The multiparameter formalism of SPDC in quantum interferometry [Ï2j has been ex- 
tended to incorporate the longitudinal nonlinearity profile of the medium. The quantum- 
interference pattern was shown to critically depend on the specific design of the nonlinearity 
profile. We studied the case of a cascaded pair of bulk crystals and experimentally verified 
the theoretical predictions. We have demonstrated that the quantum interference is sen- 
sitive to the medium between the crystals, as well as the design of the optical system for 
the down-converted light. In particular, collective interference effects were seen between the 
two probability amplitudes corresponding to detection of a photon pair generated in either 
crystal. The visibility of this collective interference between the two amplitudes depends 
on a relative phase, which is a function of the dispersion in the linear medium and the 
acceptance angle of the optical system. In the small-aperture approximation, we can con- 
tinuously sweep the quantum-interference pattern from a high-visibility triangular dip to a 
high-visibility peak, simply by changing the distance between the two crystals. In principle, 
the same effect can also be observed when the dispersion properties of the linear medium 
are changed for a fixed distance between the crystals. Furthermore, as we increased the 
aperture diameter, and thus admitted a greater range of transverse wave vectors into the 
optical system, we observed a contraction in the oscillation period of the visibility. 

Our findings are expected to be of interest to the development of SPDC sources using 
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multiple-crystal configurat ions and/or periodically poled materials, and to the advancement 
of quantum technologies through quantum-state engineering. 
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APPENDIX A: THE VISIBILITY FUNCTION V(r) 

The purpose of this Appendix is to derive the visibility function V(r) defined in Eq. fl45[) . 
We consider the special case of degenerate collinear Type-II SPDC from two cascaded crys- 
tals of identical materials but different thicknesses. This source of SPDC is used in the 
interferometer shown in Fig. 4 with symmetric apertures and no spectral filters. Given this 
explicit configuration, we evaluate Eqs. fl37j), fl38|), (|4T"|), and fl4l2| ) and rearrange the results 



to obtain Eq. (||). 

To calculate J-^b i n Eq- (f£0D, we need to consider the explicit form of the transfer function 
H for a given optical system. In the Fresnel approximation (which is well satisfied under 
the conditions of our experiments) the transfer functions Hi (i = A, B) in Eq. ([ÏO]) are given 
by 



fl-,(xi,q;a;) oc f(u) e"'^"^ 1 ) e -i&M a 

x / dy Pi (y) e <^-^ y (Al) 



where f(u) is a (typically Gaussian) spectral filter profile, u = + v is the frequency of 
the degenerate SPDC, d\ is the distance from the output plane of the second crystal to the 
aperture, d 2 is the distance between the aperture and each detector, / is the focal length 
of both lenses, and Pi(y) is the pupil function for aperture i = (A, B). In the absence of 
interference filters [f{uj) = 1], the quasi-monochromatic field approximation [y ou®) allows 
the ^-"ab functions in Eq. (HO) to become frequency- independent, as shown in Eqs. (fB3|) and 
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The wave vector mismatch functions in Eq. (|32|) and (^) are given under these approx- 
imations by 

A (q, v) = -vD + ^£ + M q (A2) 
fcp 

where A is the wave vector mismatch function due to dispersion in the crystals, and D and 
M are material properties of the crystals |22] , [23|] . Here D = u^ 1 — u' 1 is the dispersion 
coefficient, where 1i e £11*6 the group velocities for ordinary and extraordinary waves at the 
central frequency to p /2 and q = (i.e. 8 e = 6 P ). In Eq. (|A2]), M = -À- ln[n e (o;p/2, p )]e o . a . is 



the spatial walkoff vector, where e o a . is a unit vector pointing in the direction of the optical 
axis. 

Meanwhile, the dominant contribution to the dispersion function for the linear medium 
in Eq. is 

A' = k p -K D - K e (A3) 

where k p = n^cuT) uj^/c, K e = n d (üü e ) u e /c and K = n d (u ) u /c are the first-order expan- 
sions of k e and k , and ri\(u) is the index of refraction in the linear medium. Evidently the 
dispersion in the linear medium between the two crystals simply contributes a difference in 
optical path length to the overall dispersion function. For a delay line of thickness l T , the 
phase functions í/v(q, v) in Eq. ([39]) are given under the Fresnel and quasi-monochromatic 
field approximations by 



|2 



2 q 

-uD T + + M T • q — K Q — K c 



lr (A4) 



where D T and M T are the dispersion coefficient and the spatial walkoff vector for the bire- 
fringent material of the delay line, respectively. If we consider the contribution e -lí?T ( q,i ^ for 
the delay line along with Equation (p2[), we obtain 

X^ 2) (q, v) e'^' u) = e - iA ' dS ^ f dz Q r (u- z) MM- z) (A5) 

J — L r 

for r = (1,2) where xi^(q, v) and X2 : (q> u ) are the nonlinearity profiles of the first and 
second crystals, respectively. We have here introduced the quantities 
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M r (q; z) = ji^rZ-Cr)-* e 2i V |q| (A6) 
Q r (u; z) = e~ iu{DrZ+Tr) (A7) 



where the parameters 



r r = T- L 2 D ■ 5 r>1 (A8) 

p r = l T + (L 2 + d) • í r , a (A9) 

C r = M T l T + M 2 L 3 _ r • 6 rA (A10) 

for r = (1, 2) and where r = — l T D T . Under these conditions the single crystal components 
of the coincidence-count rate are then given by 

R (r) = MAo.Bg + MBo,Ae j ^ ^ ^ J[— (M r Z + £ r ), 0, S r , S r ] (All) 

and 

RÍl = 2 ^ Ao ' Bc ^ Bo ' Ac |^n Lr (z)n Lr ^-z) J[-(M r z + £ r ),Z r , Sr , Sr ] (A12) 
where for brevity we introduce U Lr (z) = rect[o,z, r ](2). Here the distance s r = d\ — z + /3 n 



Z r = -2[M r ^ + £ r ), (A13) 



and 



J[z (z),Z(z,/),s fc ,s n ] = J dqd^M k {q:,z)MM;z') T AB {q, ±c{ , u) . (AU) 

We have assumed, in keeping with our experiments, apertures which are symmetric such 
that |í»A,B(y)| = \pa,b(— y)|- Further details of .Fab(çI, 27 ) can be found in Appendix 0. 
In experimental practice, the distance d\ between the output plane of the second crystal 
and the aperture is much longer than the crystal thicknesses L\ and L 2 and the optical path 
length l T of the delay line. As such, we may assume that 

s r = dt - z + p r ~ di + (2 - r)d. (A15) 
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Also note that there are no approximations constraining the value of d. 

The baseline coincidence-count rate from collective interference between the amplitudes 
from the two crystals is given by 



2 2 
D (12) /^Ao.Be + /%o,Ae 

=e 



■ \ I dz n L2 [z-IL—HJ X[-(M 1 z + C 1 ),2° 12 (z),s 1 ,s 2 } e~ iA ' d + 

+ I dz Ii/.j :) 11^ ( z - Z[-(M 2 2 + £ 2 ), 2>\\{z), s 2 , si] e iA ' d 



(A16) 



where 



^(z) 



2 (Mi - M 2 ) + M 2 ^^ + d - C 2 



D 



(A17) 



and the remaining parameters are given by 



n = r - L 2 L> si = di + d Ci = M T l T + M 2 L 2 
t 2 = t s 2 = d\ C 2 = M r / T . 



(A18) 
(A19) 



It is important to note that Rq = for all vàlues of r, since in this case T\ — r 2 = —LD 2 
and the two rectangular functions Hi, r {z) in Eq. ( |A16| ) never actually overlap. Thus the 
coincidence-count rates are strictly constant outside the region of interference. This lack of 
shoulder modulation is an important indicator of the purity of the observed polarization- 
entangled two-photon state. The collective interference term itself is given by 



R (12) _ 9 H-Ao,Be ^Bo,Ae 



dzU Ll (z) U L2 [H±H- Z \ l[-(M 1 z + C 1 ),Z^(z),s 1 ,s 2 ] e-^' d + 



dz IÍ L2 (z) Ii Ll 
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íA'd 
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z \ I[-(M 2 z + £ 2 ),Zff{z),s 2 , Sl ]e 



(A20) 



with 



z (M x - M 2 ) + M 2 ^-±p + Cx + C 2 



where Z 2 Í nt can be found by interchanging all indices (1 <-> 2) in Z 



0,int 
12 • 



(A21) 
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With sufíicient àlgebra, it is possible to arrange all of the contributions to the coincidence 
rate in Eqs. (|MID , (|XI2|) , ( |M6| ) and ( |A20l ) to obtain the structure of Eq. (H), where the 
baseline coincidence-count rate is given by 



2 r 



*= T. 



L\ L2 



Pa(0) Pn(0) 



Dsj Ds 2 2 

where -Pi(q) is given by Eq. (|B4] ), and the polarization analyzer projection factor is 

r, /^Ao,Be /^Bo,Ae 



(A22) 



^pol 



^Ao.Be + ^BoM 



(A23) 



This, in turn, leads to Eq. 
L\ = L2 = L: 



the desired visibility function in the special case where 



+ 1 / (fa n lW n 1 (|- ï ) 



z r 



-2e 



/2r 



3 / ^n L (^)n L ^— -l-z) h\ 



eM(£ + i ; JL 



. (A24) 



Here the functions (? r= i 2 are given by 



Z T 



G (r) ( ^ ) = [-(M r z + C r ),Z r , s r , s r 



(A25) 



for the single-crystal contributions, while the component due to collective interference is 
given by 



{- + -; — 

L 2 LD 



■ Àí 



-(M l z + Ci),Z° 12 (z),si,s 2 



(A26) 



The function ÀÍ is defined in Appendix [B] by Eq. ( p9|) , while the normalization factor is 
given by p = (d\ + d)/d\. 

For a delay line comprised of a thin z-cut birefringent element such as quartz, M T = 
and 



M 2 L 2 , 



0. 



(A27) 



Under this assumption there are two important limits: large separation between the crystals 
(d — > 00) and contact between the crystals (d — > 0). When the two crystals are moved very 
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far apart, we expect interference from SPDC to be governed only by the second crystal, and 
the quantum-interference pattern to be identical to that of SPDC from a single crystal of 
thickness L2. This can be seen from Eq. (|A24|) by noting that in this case the normalization 
factor p —>■ 00 and only the last of the three integrals in Eq. ( |A24| ) survives. In this limit 

^oo(r) = jjdz U L (z) n L - *) (£; JL) (A28) 

and the shoulder normalization factor ( |A22|) is identical to that of a single crystal of thickness 
L. When two crystals of the same material with parallel optical axes are in contact, the 
quantum interference pattern is identical to that obtained from a single crystal of thickness 
2L. However, note that this does not hold for antiparallel orientations of the optical axes of 
the two crystals. 



1. Parallel Optical Axes 

If the crystals are oriented such that their optical axes are parallel, Mi = M 2 = M and 
the Q functions in Eq. (f45|) are specified by 

x) = Àí [-ML(C + x), -2MLx, d x + d, d x + d] (A29) 
gl 12) (C;x) = AÍ[-ML(Ç + x),-2MLx,d 1 + d,d l ] (A30) 
gW((;x)=N[-ML(Ç + x),-2MLx,d 1 ,di] (A31) 

where ( = z/L, x = r/LD, and Aí is defined in Eq. (p9| ). 



2. Antiparallel Optical Axes 

The antiparallel orientation of the optical axes of the two crystals gives M x = — M 2 = M 
and the ^-functions are specified by 

gW(( ; X )=ÀÍ [-ML(C + x - 2), -2ML(x -l),d 1 + d,di + d] (A32) 
g£ 2 \(; x)=ÀÍ [-ML(C + x - 2), -2ML(C - 1), d x + d, tZj (A33) 
gW{(;x) =ÀÍ[ML{C + x),2MLx,d 1 ,d 1 ]. (A34) 
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APPENDIX B: THE FUNCTION N 



In this Appendix, we derive an explícit form of the function ÀÍ, which appears in 
Eq. ( \K2Üj ). We define 

J(z ,Z, 



MM H , Sn ) (V2)2 pMo) 

where Pj(q), k,n— (1,2), si and s 2 are defined in Eqs. ( |A18| ), and 



(Bl) 



l[z {z),Z{z,z'),s k ,s n } = J d<idc[M k {^z)M* n {o[]z') .F AB (q,±q», (B2) 

where M. r is defined in Eq. (|A6|). The function JF AB , defined by Eq. (PXf), under the Fresnel 
approximation and in the absence of interference filters takes the form 



•F Ais(q, ±q', u) = P A (q T q') ft(-q ± q') e" 1 ^^ 1 Hq 1 > . 



(B3) 



Here 



A(o)= / ^y^(yM(y)e- iy - q 



di(f) = dx 



for i = (A, B). Given these explicit forms, Eq. ( |B2| ) becomes 



X(z , Z, Sjt, s r , 



(fcp/2) s 



exp 



-i-^|Z| 2 

8Sn 



(B4) 
(B5) 



dq P A (-q) P B (q) W í q+ Zj e"^ 20 

(B6) 



where 



W(q) 



2idí fcn) 



exp 



2idi fcn ) 



q 



with 



(B7) 



d! 



(fen) 



Sfe 



Finally, using Eq. ( p6|) we obtain 



(B8) 
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-i^|Z| 2 

8s n 



dq. P A (-q) P B (q) W í q + ^Zj e" 1 ^ (B9) 



A/"(z , Z, s fe , s n ) = exp 
with Vi(q) = -Pj(q)/-Pj(0) for i — A, B. In the limit s fc — > s n , Eq. ( |BÜD becomes 
A/"(z , Z, s n , s n ) = exp 



-i-^|Z| 2 

Bs n 



4s n 



(B10) 



1. Small-Aperture Approximation 

Suppose that the apertures are so small that the pupil functions í>A,B(y) ma Y be treated 
as delta functions at y = 0. In this case, Va(<í) is a constant and Eq. ( |B9|) becomes 



A/" SM (zo, Z, s n ) = exp 
In the limit Sk — > s n we then have 

JV S m(zo,Z, s n ,s n ) = exp 



— í- 



.fcn /|Z - Z | 2 |Z | 2 ' 



(Bll) 



-i^íiz 

8s n V 



|2 I |2 

zo - z 



(B12) 



2. Gaussian-Aperture Approximation 

To simplify the analysis, we consider the case of Gaussian apertures with pi{y) 
exp(— |y| 2 /2r 2 ) and Vi(q) = exp(— |q| 2 rj/4), so that 



A/g(zo, Z, s k , s r 



1 — Z7 

1 , 2 eX P 

1 4- y 
x exp 



A't. 



— 1 



8(s fe - s n ) 



zo — Z 



where 



7 " U kn) 2 
It is useful to write the complex constant 

1 — ry 1 

1 + 7 2 ~ V^+Y ' 
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-i arctan 7 



(B13) 



(B14) 



(B15) 



so that we may write the phase 



7 = — arctan 



K r l + r l d/ di 



Adx 2 l + d/dj' 
In summary, the A/"-function in the Gaussian-aperture approximation is given by 



■A/g(z , Z, s k , s n ) 



1 



VTT7 



f ex p 



À r- 



7 



x exp 
x exp 

In the limit Sk — > s n , we have 

J\í G (z ,Z,s n ,s n ) = exp 



À' T . 



8d (kn) 1 + 7 2 
1 



— 1 



8d (kn) J + 7 2 



z - 



Sk 



Sk S n 
2" 

-z 



8(Sfc - Sn) 



|Z| 2 + i(^ 



A; P Z 



4s r , 



2 2 i 2 



x exp 



-^(iZ-^-lzol 2 ) 



(B16) 



(B17) 



(B18) 
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FIGURES 

FIG. 1. Control over the nonlinearity profile of the gener ation medium allows control over the 
SPDC state function <J?(A). In the case of a monochromatic plane wave pump, $(A) is simply 
the inverse Fourier transform of the nonlinearity profile x^ 2 \ z )- The upper figure shows a sketch 
of |<Ï>(A)| 2 for Type-II SPDC from a single bulk crystal of thickness L. The lower figure shows a 
sketch of |$(A)| 2 for SPDC from a crystal with sinusoidally varying nonlinearity. In principle, any 
weighting profile of the signal and idler photons can be obtained by a judicious choice of crystal 
structure. 

FIG. 2. Impact of generation-medium symmetry on the weighting profile of SPDC pumped by 
a monochromatic plane wave. The upper figure shows |<3?(A)| 2 for SPDC from two bulk crystals 
of unequal thickness separated by dispersive linear medium. The lower figure shows |$(A)| 2 for 
SPDC generated by a cascade of two bulk crystals of equal thickness. 

FIG. 3. Propagation of the signal and idler photons through arbitrary linear optical systems 
TÍAj and W Bj • 

FIG. 4. Schematic of a polarization interferometer for which we compute quantum-interference 
patterns. In this illustration, collinear SPDC is generated in two bulk crystals of arbitrary thickness 
separated by an air gap. 

FIG. 5. Sketch illustrating a heurístic approach to calculating quantum interference in the 
single-mode limit for two bulk crystals of the same material and thickness separated by a linear di- 
electric. The results are coincidence rates as shown in the bottom right inset. The coincidence-rate 
data in the upper left inset is from representative experiments. Details are provided in the text. 

FIG. 6. A figure analogous to that of Fig. 5 for the case of two bulk crystals of unequal thickness. 
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FIG. 7. Visibility |V(LD)| as a function of crystal separation d for the case of parallel optical 
axes. The curves are plots of Eq. ( |54| ) for aperture diameters b = 2.5 mm, 4.0 mm, and 5.0 mm. 
Sòlid, shaded, and open circles denote the màxima of visibility. The distance d\ (see Fig. 4) is 750 
mm. Note that the period of visibility modulation as a function of crystal separation d contracts 
for increasing aperture diameter. 

FIG. 8. Visibility at r = LD for two identical crystals in contact (d = 0) as a function of the 
(1/e) width r of identical Gaussian apertures in each arm of the interferometer shown in Fig. 4. 
Sòlid (dotted) curves correspond to parallel (antiparallel) optical axes. Results for different crystal 
thicknesses L are shown concurrently. The curves were generated with the parameters d\ = 1 m, 
A p = 2ir/k p = 351.1 nm and |M| = 0.07. 

FIG. 9. Experimental apparatus used to study polarization interference in collinear SPDC from 
two cascaded crystals separated by an air gap of thickness d. Details are found in the text. 

FIG. 10. Visibility function at the center of the interference region, V(LD), as a function of 
crystal separation d. Symbols represent data from parallel (shaded circles) and antiparallel (open 
circles) optical axes. Note that in our experiments V(LD) = 1 corresponds to an interference 
dip, while V(LD) = — 1 corresponds to an interference peak. Insets at top are representative 
interference patterns at extremes of V(LD), taken with d = 17.5 mm and d = 32.5 mm. Sòlid 
curves are plots of Eqs. (54) and (|6^) for parallel and antiparallel orientations, respectively. 



FIG. 11. Visibility function V{LD) as a function of crystal separation d for three different 
aperture diameters in the case of parallel optical axes. Sòlid curves are plots of Eq. (|54|). Sòlid, 
shaded, and open circles are experimental data for aperture diameters b = 2.5 mm, 4.0 mm, and 
5.0 mm, respectively. The distance d\ (see Fig. 4) is 750 mm. 
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